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Abstract 

We study 4a linear-chain states in 16 O in comparison with 3a states in 12 C by 
using the Generator Coordinate Method within a microscopic ./Va-cluster model. 
It is shown to be very important to solve relative motion of a-clusters by taking 
into account the orthogonality between the linear-chain states and other low- 
lying states including the ground state. For 12 C, the 3a chain state disappears 
because of this orthogonality. The 4a chain state in 16 0, on the other hand, is 
hardly affected by low- lying states and persists to remain above the 4a threshold. 
The calculated moment of inertia of the 4a linear-chain rotational band (88 keV) 
reproduces the experimentally suggested value qualitatively. 



* E-mail: itagaki@nucl.phys.hokudai.ac.jp, Fax: +81-11-746-5444. 



1 Introduction 



A nuclear study of a large deformation is one of the frontier in nuclear spectroscopy. 
Since the first observation of the superdeformed band |TJ , various superdeformed nuclei 
(~ 1:2 deformation) have been found, and nowadays extensive studies on these bands 
including their excited bands are in progress. Furthermore, there are several reports 
on the hyperdeformed band in which the deformation is around 1:3 0. These large 
deformed nuclei are mainly found in the medium-heavy mass number region where 
the liquid-drop energy and the shell-correction balances. Therefore, larger deformation 
than 1:3 is not expected to appear in this mass number region, since the shell-correction 
is not enough to make a barrier at larger deformations. 

On the other hand, a large deformation above 1:3 may be possible in light nuclei 
due to the a-cluster structure. For example, the first example of the "superdeformed" 
nuclei is nothing but 8 Be which has the a — a structure, and there have been a lot 
of discussion on the nature of state of 12 C as the 3a linear-chain state ||, |J. Ex- 
perimentally, it is more difficult to observe large-deformed states or linear-chain states 
of a-clusters in light nuclei, since these states have large particle decay widths, and 
consequently a coincident detection of multi-a particles becomes necessary. From a the- 
oretical point of view, on the other hand, it is required to estimate the couple-channel 
effects over the broad range of deformation or geometric structure, since a mean field 
picture may not be valid for the highly excited states in light nuclei. However, since it 
gives us the opportunity to explore the most deformed nuclei, the study of these Na 
linear-chain states || is valuable. 

For example, there are several experimental candidates in light 4N nuclei for 
large deformed states above 1:3. The first example is the 4a linear-chain band starting 
around the 4a threshold energy region in 16 suggested by Chevallier et al. || through 
the data analysis of the 12 C+a — > 8 Be+ 8 Be reaction. This suggestion is supported by 
a theoretical work by Suzuki et al. 0; they analyzed the decay widths from such a 
strongly coupled linear-chain state, and discussed that the excited states observed by 
Chevallier et al. are possibly characterized by a 4a linear-chain structure. The second 
example is the 6a linear-chain states, which is discussed as a candidate to explain a new 
molecular resonance state observed by Wuosmma et al. M in the excitation function 
of the inelastic reaction 12 C( 12 C, 12 C(0^)) 12 C(0^) at E X =5GA MeV 0. Although the 
clarification of the existence of a 6a linear-chain state is very challenging, the coupling 
to 12 C*- 12 C* and 16 0*- 8 Be may be important QTOj and makes it difficult to study a 
pure linear-chain state. Therefore, it is important and urgent to study the stability of 
linear-chain states in 12 C and 16 O at this stage. 

The reliability of the linear-chain structure has been usually studied through 
kinematical analyses of decay widths. It has been also discussed on the analyses of 
small vibration around the equilibrium configuration by Ikeda |TTJ. In those previous 



studies, however, the linear-chain structure was assumed a priori. Namely, the coupling 
effects with other states have not been taken into account. These couplings, especially 
with lower excited states are expected to play an important role on the stability of 
linear-chain states, since they push-up the energy of linear-chain states and generally 
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increase the instability. Therefore, it is necessary to calculate in a wide wave function 
space, which covers not only the linear-chain configurations but also the lower excited 
states. 

In this paper, we investigate the stability of 4a linear-chain configurations in 
16 comparing with 3a states of 12 C by applying a microscopic iVa-cluster model. 
In contrast to the previous works, we discuss the stability of the Not linear-chain 
structure by solving dynamics of Not systems. This is achieved by diagonalizing the 
Hamiltonian matrix in a wide space which also covers low-lying levels including the 
ground state. The framework used here is based on the Generator Coordinate Method 
(GCM) by which we can solve relative motion of a clusters by taking account of the 
superposition of many different intrinsic configurations. We use the Constraint Cooling 
Method jl2 , 13 , 14, [15| proposed in the framework of AMD [16|] in order to generate 



GCM basis which describes low-lying levels. This method has been already shown to 
provide the suitable GCM basis states effectively p 



This paper is organized as follows; in section 0, we explain our model especially 
the way to construct the GCM basis states. In section ||, we show our results that the 
12 C linear-chain state disappears because of the orthogonality to other states, however 
the 16 linear-chain state exists above the 4a threshold and shows the rotational band 
structure. In section |] summary and conclusion are presented. 



2 Method 

In the framework of GCM, the energy spectra and the wave functions are solved by 
diagonalizing the Hamiltonian matrix through the Hill- Wheeler equation; 

V I //U X'-, <Ir / ' / • (!) 

j 

J2((^ M \H\H M ) - E(^ M \H M })c k = , (2) 

k 

where <f>J M is the jth GCM basis wave functions with the total angular momentum J 
and its z-component M. Since the reliability of the GCM calculation largely depends 
on the adopted basis wave functions, we have adopted the GCM basis states which 
reasonably reproduce the low-lying excited states in addition to the linear-chain con- 
figurations described in the next section. The first set of GCM basis states is generated 



by applying the constraint cooling method |13] to the Bloch-Brink wave function for an 



Na system. The angular momentum projection from these wave packets is carried out 



numerically. A detailed explanation of this method is found in our previous paper [17 



Hereafter, we briefly summarize our method. The Bloch-Brink wave function for 



an Na system is the Slater determinant of Gaussian wave packets [12, 16 1 with the 
assumption of the [4]-symmetry for the spin and iso-spin part, 



A' 



$(K;Z l }) = ^nilC(rr;Z l ), (3) 



i=l a=l 
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ft (r Q ; Z«) = exp{-u(r a - Z,^) 2 + Z 2 /2} Xa (i = 1 ~ iV), (4) 

where x a = {pf, p|, nf, n|} are the spin-isospin wave functions. In the spatial part, 
the complex parameters Z$ = v / z/(Dj + iK.i/2vK) represent phase space coordinates. 
In this work, the size parameter v = 0.27 fm -2 is chosen to reproduce the energy and 
the r.m.s. radius of the a particle. 

In order to generate various intrinsic states, we employ the constraint cooling 
method fll2] , |TB| , II, |T3|, by which we can search for a set of suitable complex parameters 



{Zi} to minimize the expectation value of the energy under some constraints. Following 
Ref. ||13[| , we obtain an optimum configuration state by solving the constraint cooling 



equation given by 

dz 1 = _m d_m 

d(5 dZ* fr[ 11 dZ* ' 1 ; 

where TC is the expectation value of the Hamiltonian operator. The expectation value 
of the constrained observable (Wi) in Eq. (|5]) plays a role of the generator coordinate, 
and t)i is the Lagrange multiplier to keep the given constraint along the evolution. In 
this work, we constrain the expectation value of the total oscillator quantum numbers 
Wi = n x , n y and n z in a rectangular coordinate system JHJ. It is possible to constrain 



other observables, such as an angular momentum, however, cooling with an angular 
momentum constraint essentially results in giving the rotating nucleus with the similar 
intrinsic state when the angular momentum is not so high. Then these states are not 
very effective as GCM basis states. 

Since the resultant cooled state is not the eigenstate of an angular momentum 
in general, it is necessary to carry out the angular momentum projection to construct 
GCM basis wave functions in Eq. ([[]). 

o J _|_ i t 

$^({rf ; Z t }) = J dQD{* K (m(m({rt; Z t }), (6) 

where Q is the Euler angle of {r"}. In the numerical integration over the Euler angle 
f2(a,/?,7), we have adopted 25 3 mesh points and integration is carried out with the 
Gauss- Legendre integral technique. The intrinsic angular momentum component K is 
chosen to be zero for simplicity. 



3 Na linear-chain states in C and b O 
3.1 3a states in 12 C 

In many previous analyses, it has been shown that a large part of low-lying excited 
states of 12 C have clustering structures. Although Morinaga proposed the linear- 
chain structure for 12 C, the clustering states of O2 (7.66 MeV) and 22~(~10 MeV) have 
been discussed to have a loosely-coupled 3a configuration rather than the linear- 
chain structure. We here investigate the stability of the 3a-linear-chain configuration 
in 12 C from a viewpoint of orthogonality to low-lying states by applying a microscopic 
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Figure 1: The linear-chain configurations of 3a (a) and 4a (b). 

3a model. In this calculation, we use the Volkov No. 2 potential Jl8| with the Majorana 
value M = 0.592 as an effective nucleon-nucleon interaction. By using this interaction, 
the experimental binding energy of 12 C (7.4 MeV) from the 3a threshold, which is 
calculated as —83.69 MeV, is well reproduced. 

Various intrinsic 3a configurations can be described by a parameter set (Zi, I2, 9), 
where l\ and I2 are intervals between two a-clusters and 9 is a bending angle (See Fig.l 
a)). The linear-chain configuration corresponds to 9 = 71. Here, we treat l\ and 9 as 
generator coordinates and construct the basis states which have various sets of these 
coordinates. We employ li = 1,2,3,4,5 fm and 9 = n ~ n/3. We regard l 2 as a 
variational parameter and apply a cooling method to search for the energy minimum. 
The 45 GCM basis states described by parameters 9) are obtained, and their 

explicit values are presented in Table I. 

After projecting the obtained pure linear-chain intrinsic states (9 = n) to J n = 
+ , we obtain the energy —71.72 MeV at the minimum point (li, I2) =(4 fm, 2.89 fm) 
(GCM basis No. 4). This energy is much higher than the 3a threshold where the 0^ 
state is observed. To solve the longitudinal vibration of three a clusters, we employ 5 
GCM basis states which are given by I2) = (1 fm, 3.55 fm), (2 fm, 3.22 fm), (3 fm, 
2.97 fm), (4 fm, 2.89 fm), and (5 fm, 2.92 fm), keeping 9 = ir (GCM basis of No.l ~ 
No. 5). After solving the GCM equation, we have the binding energy gain by 2.86 MeV 
(—74.59 MeV), but its energy is still too short for the energy of the experimentary 
observed state. 

Furthermore, we take account of the bending motion (9 = ll7r/12, 107r/12 and 
97r/12) which corresponds to fluctuation around the linear-chain configuration. The 
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Table 1: The GCM basis states for 12 C described by Z x (fm), l 2 (fm) and 9 (radian) 
of Fig. 1 (a), li is determined to minimize the energy at given l\ and 9. 

solutions of the GCM equation with those basis states (No.l ~ No.2U) give the lowest 
energy — 8U.4U MeV. Although this energy is about 3 MeV higher than the energy of 
the 0% state, it seems to give a good correspondence with that. The wave function of 
this solution has a dominant components of the linear-chain configuration, which are 
seen from the fact that overlap with the linear-chain GCM solution (described in the 
previous paragraph) is U.65. 

From the analyses described above, one may consider that the 0% state can be 
interpreted as a 3a linear-chain state with small bending vibrations. However, this 
solution may also have a large overlap with the ground state of the 3a system, and it is 
necessary to solve the U2" state by taking into account of the coupling with the ground 
state. In order to clarify this point, we performed the GCM calculation with larger 
basis states including those which describe low-lying levels. After diagonalizing the 
Hamiltonian with GCM basis states (No.l ~ No. 45 in table I), we obtain the ground 
state energy —91.23 MeV which well corresponds to the experimental one (—92.2 MeV). 

In Fig. 2, we show the conversion behavior of the energy spectrum. As the number 
of GCM basis states increases, each basis configuration generally becomes closer to 
equilateral triangle states. The ground U + state is reproduced around the basis states 
No.3U which corresponds to the configurations of 9 = 6tt/12. By this increase of 
basis states, the second U + state converges energetically. This solution has the overlap 
with linear-chain {9 = tt) basis states of only U.346 (absolute square), and cannot be 
interpreted as a linear-chain state. 
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Figure 2: The diagonalization of Hamiltonian of 12 C. 



3.2 4a states in 16 



In a similar way to the 3a system, we investigate the stability of the 4a linear-chain 
configuration in 16 0. We use the Volkov No. 2 interaction [|18| with the Majorana pa- 
rameter M = 0.63, in order to reproduce the ground state energy of 16 (— 127.6MeV) 



17]. However, this effective nuclear interaction has a shortcoming that the calculated 
binding energy of 12 C (-82.0 MeV) does not reproduce the observed one (-92.2 MeV). 
This problem has been discussed jL7j from both sides of model spaces and effective in- 
teractions. As one important conclusion, it has been shown that the clustering states 
are well reproduced around the corresponding threshold energy for various effective 
interactions. 

We apply the microscopic 4a-cluster model and carry out the GCM calculation 
for stability examination of the 4a linear-chain configuration. The employed GCM 
basis states are listed in Table II. The linear-chain states and their bending motions 
are described by 45 GCM basis states (No.l ~ No.45 in Table II) of the 8 Be+ 8 Be 
configurations with length parameters (l\ and Z 2 ) and an angle parameter (8 = 71 ~ tt/2) 
as shown in Fig.l (b). The relative position of two 8 Be clusters is determined by / 2 
(Z 2 = 0.5, 1.0, 1.5,2.0 and 2.5 fm) and the distance (Zi) between the two a-clusters of 
each 8 Be is chosen to minimize the energy at given Z 2 and 9. In addition, we prepare 
15 GCM basis states obtained by the Constraint Cooling Method for the description 
of other lower states. 

As the first step in research of the 4a-linear-chain state, we search the + optimal 
linear-chain configuration with (li = 3.29 fm, l 2 = 1.50 fm, 6 = 0) at the energy 
E = 34.79 MeV from the 4a threshold within a single Slater determinant. As the next 
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Table 2: The GCM basis states for 16 0. Basis state configurations configuration No.l 
~ No. 45 are described by l\ (fm) k (fm) and 9 (radian) of Fig. 1 (b). The value of I2 
is determined to minimize the total energy. The GCM basis states No. 46 ~ N0.6U are 
obtained by the constraint cooling method, and classified with the constraint principal 
quantum numbers (n x , n y and n z ). 



step, we examine its stability for mixing with configurations other than the linear-chain 
by solving the GCM equation with increasing the number of the basis states. 

The U + -state energy convergence is shown in Fig.3, and the ground state energy 
converges to —129.1 MeV. The energy of the U + state obtained at —91.68 MeV (finally 
corresponds to U^j~) converges with a very small number of GCM basis states (= No. 15), 
which means this state contains large linear-chain components with small bending 
(< IUtt/12). As an evidence, this state has an overlap (absolute square) with the GCM 
basis states Nos.l-lU (9 = ix and Htt/12) of U.461 and with the GCM basis states 
Nos.1-15 (6 = ix, Htt/12 and IUtt/12) of U.613. Even when a larger number of GCM 
basis states are adopted, the state of —91.68 MeV is hardly affected by other lower 
states and approximately keeps its own energy value. However this energy is higher 
about 2U MeV compared with the 4a threshold energy (—111.6 MeV). 

There may be at least two reasons for the too-high-energy solution of U^ : 1) The 
12 C+a threshold is not correctly reproduced due to the present iViV interaction, and 
is calculated at the energy lower than experiment by about 1U MeV. As was discussed 
in the previous paper, it is very difficult problem to reproduce both threshold energies 
of 12 C+a and 4a by an effective A^A^ interaction. Because of this shortcoming, 4a 
clusterized configurations are calculated at rather high excited energies. 2) The GCM 
basis states describing the Ug may be insufficient. We have employed the configurations 
of bending motion deviating from the linear-chain, and of the three-dimensional 4a 
structure though they contribute little to this Ug" state. Then the energy of this state 
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Figure 3: The diagonalization of Hamiltonian of 16 0. 



would be lower if we solve other small vibrations around the linear-chain configuration 
such as Z letter-type mode or three-dimensional twinst mode. 

However, an improved calculation based on the above reason for the fault would 
only increase the quantative precision, and at least we can suggest two things from the 
resent result. First, the state of —91.68 MeV has a large overlap with the linear-chain 
configuration which shows the 6 stability. Second, this state is hardly affected by other 
lower energy states and there is a large possibility to correspond to remain as resonance 
state. 

In Fig. 4, the linear-chain states with different angular momenta are presented. 
From the 2 + state, the rotational band is shown to appear with a very small slope of 
88 keV because of its large inertia moment. This moment of inertia is similar to the 
experimentally suggested value (about 60 keV). 

4 Summary and Conclusion 

In this peper, we have microscopically studied the existence and the stability of Na 
linear-chain state in 12 C and 16 0. We have emphasized the importance of the orthog- 
onality between the linear-chain state and other lower states. The linear-chain state 
of 12 C is hardly considered to be a stable state as we increase the number of the 
GCM basis states. In other words, their components fragment to the levels over broad 
range of excitation energy, especially to the 0^~ state, and these couplings make a pure 
linear-chain state unstable. On the other hand, the linear-chain state of 16 is hardly 
affected by other states, and the rotational band structure is shown to appear, although 
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Figure 4: The rotational band of 16 linear-chain state. The numbers show the order 
of excited state. 



the absolute energy of the band head is higher by about 20 MeV compared with the 
4a-threshold. 



These conclusions are consistent with the previous work by Ikeda ||TT|, where the 
stability against small vibrations around the equilibrium configuration is considered. 
However, we have shown in this work that the linear-chain state of 16 may survive 
even after the couplings to the lower states are taken into account. 
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Abstract 

We study 4a linear-chain states in 16 in comparison with 3a states in 12 C by 
using the Generator Coordinate Method within a microscopic iVa-cluster model. 
It is shown to be very important to solve relative motion of a-clusters by taking 
into account the orthogonality between the linear-chain states and other low- 
lying states including the ground state. For 12 C, the 3a chain state disappears 
because of this orthogonality. The 4a chain state in 16 0, on the other hand, is 
hardly affected by low-lying states and persists to remain above the 4a threshold. 
The calculated moment of inertia of the 4a linear-chain rotational band (88 keV) 
reproduces the experimentally suggested value qualitatively. 
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1 Introduction 



A nuclear study of a large deformation is one of the frontier in nuclear spectroscopy. 
Since the first observation of the super deformed band [1], various superdeformed nuclei 
(~ 1:2 deformation) have been found, and nowadays extensive studies on these bands 
including their excited bands are in progress. Furthermore, there are several reports 
on the hyperdeformed band in which the deformation is around 1:3 [2]. These large 
deformed nuclei are mainly found in the medium-heavy mass number region where 
the liquid-drop energy and the shell-correction balances. Therefore, larger deformation 
than 1:3 is not expected to appear in this mass number region, since the shell-correction 
is not enough to make a barrier at larger deformations. 

On the other hand, a large deformation above 1:3 may be possible in light nuclei 
due to the a-cluster structure. For example, the first example of the "superdeformed" 
nuclei is nothing but 8 Be which has the a — a structure, and there have been a lot 
of discussion on the nature of Oj" state of 12 C as the 3a linear-chain state [3, 4]. Ex- 
perimentally, it is more difficult to observe large- deformed states or linear-chain states 
of a-clusters in light nuclei, since these states have large particle decay widths, and 
consequently a coincident detection of mult i- a particles becomes necessary. From a the- 
oretical point of view, on the other hand, it is required to estimate the couple- channel 
effects over the broad range of deformation or geometric structure, since a mean field 
picture may not be valid for the highly excited states in light nuclei. However, since it 
gives us the opportunity to explore the most deformed nuclei, the study of these Net 
linear-chain states [5] is valuable. 

For example, there are several experimental candidates in light AN nuclei for 
large deformed states above 1:3. The first example is the 4a linear-chain band starting 
around the 4a threshold energy region in 16 suggested by Chevallier et al. [6] through 
the data analysis of the 12 C+a — > 8 Be+ 8 Be reaction. This suggestion is supported by 
a theoretical work by Suzuki et al. [7]; they analyzed the decay widths from such a 
strongly coupled linear-chain state, and discussed that the excited states observed by 
Chevallier et al. are possibly characterized by a 4a linear-chain structure. The second 
example is the 6a linear-chain states, which is discussed as a candidate to explain a new 
molecular resonance state observed by Wuosmma et al. [8] in the excitation function 
of the inelastic reaction 12 C( 12 C, 12 C(0+)) 12 C(0+) at E X =56A MeV [9]. Although the 
clarification of the existence of a 6a linear-chain state is very challenging, the coupling 
to 12 C*- 12 C* and 16 0*- 8 Be may be important [10] and makes it difficult to study a 
pure linear-chain state. Therefore, it is important and urgent to study the stability of 
linear-chain states in 12 C and 16 at this stage. 

The reliability of the linear-chain structure has been usually studied through 
kinematical analyses of decay widths. It has been also discussed on the analyses of 
small vibration around the equilibrium configuration by Ikeda [11]. In those previous 
studies, however, the linear-chain structure was assumed a priori. Namely, the coupling 
effects with other states have not been taken into account. These couplings, especially 
with lower excited states are expected to play an important role on the stability of 
linear-chain states, since they push-up the energy of linear-chain states and generally 
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increase the instability. Therefore, it is necessary to calculate in a wide wave function 
space, which covers not only the linear-chain configurations but also the lower excited 
states. 

In this paper, we investigate the stability of 4a linear-chain configurations in 
16 comparing with 3a states of 12 C by applying a microscopic iVa-cluster model. 
In contrast to the previous works, we discuss the stability of the Net linear-chain 
structure by solving dynamics of Net systems. This is achieved by diagonalizing the 
Hamiltonian matrix in a wide space which also covers low-lying levels including the 
ground state. The framework used here is based on the Generator Coordinate Method 
(GCM) by which we can solve relative motion of a clusters by taking account of the 
superposition of many different intrinsic configurations. We use the Constraint Cooling 
Method [12, 13, 14, 15] proposed in the framework of AMD [16] in order to generate 
GCM basis which describes low-lying levels. This method has been already shown to 
provide the suitable GCM basis states effectively [17]. 

This paper is organized as follows; in section 2, we explain our model especially 
the way to construct the GCM basis states. In section 3, we show our results that the 
12 C linear-chain state disappears because of the orthogonality to other states, however 
the 16 linear-chain state exists above the 4a threshold and shows the rotational band 
structure. In section 4 summary and conclusion are presented. 



2 Method 



In the framework of GCM, the energy spectra and the wave functions are solved by 
diagonalizing the Hamiltonian matrix through the Hill- Wheeler equation; 

y\r ni ^> <l>; /w . (1) 

3 

) - E(^ M \^ J k M ))c k = , (2) 

k 

where is the jth GCM basis wave functions with the total angular momentum J 
and its ^-component M. Since the reliability of the GCM calculation largely depends 
on the adopted basis wave functions, we have adopted the GCM basis states which 
reasonably reproduce the low-lying excited states in addition to the linear-chain con- 
figurations described in the next section. The first set of GCM basis states is generated 
by applying the constraint cooling method [13] to the Bloch-Brink wave function for an 
Net system. The angular momentum projection from these wave packets is carried out 
numerically. A detailed explanation of this method is found in our previous paper [17]. 

Hereafter, we briefly summarize our method. The Bloch-Brink wave function for 
an Net system is the Slater determinant of Gaussian wave packets [12, 16] with the 
assumption of the [4]-symmetry for the spin and iso-spin part, 

N 4 

$(K;Z 8 }) = ^nn^(rr;Z 8 ), (3) 

i=l a = l 
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<f>f(r a ; Z,-) = exp{-u(r a - Z % j^vf + Z 2 /2} X „ (i = 1 ~ N), (4) 

where \a = {p T? P I? n T? n 1} are the spin-isospin wave functions. In the spatial part, 
the complex parameters Z 4 - = y/v(Di + zK 8 /2z//i) represent phase space coordinates. 
In this work, the size parameter v = 0.27 fm~ 2 is chosen to reproduce the energy and 
the r.m.s. radius of the a particle. 

In order to generate various intrinsic states, we employ the constraint cooling 
method [12, 13, 14, 15], by which we can search for a set of suitable complex parameters 
{Z;} to minimize the expectation value of the energy under some constraints. Following 
Ref. [13], we obtain an optimum configuration state by solving the constraint cooling 
equation given by 

dz 1 = _m d^, 

d/3 dZ* ^ " 8Z* ' 1 ; 

where H is the expectation value of the Hamiltonian operator. The expectation value 
of the constrained observable (Wi) in Eq. (5) plays a role of the generator coordinate, 
and j]i is the Lagrange multiplier to keep the given constraint along the evolution. In 
this work, we constrain the expectation value of the total oscillator quantum numbers 
W\ = n X} n y and n z in a rectangular coordinate system [17]. It is possible to constrain 
other observables, such as an angular momentum, however, cooling with an angular 
momentum constraint essentially results in giving the rotating nucleus with the similar 
intrinsic state when the angular momentum is not so high. Then these states are not 
very effective as GCM basis states. 

Since the resultant cooled state is not the eigenstate of an angular momentum 
in general, it is necessary to carry out the angular momentum projection to construct 
GCM basis wave functions in Eq. (1). 

? J _|_ 1 f 

$ JM (K ; z t }) = J d nD^n)R(m({rf-, z,-}), (6) 

where O is the Euler angle of {r"}. In the numerical integration over the Euler angle 
Q(a, f3, 7), we have adopted 25 3 mesh points and integration is carried out with the 
Gauss-Legendre integral technique. The intrinsic angular momentum component K is 
chosen to be zero for simplicity. 

3 Na linear-chain states in 12 C and ie O 

3.1 3a states in 12 C 

In many previous analyses, it has been shown that a large part of low-lying excited 
states of 12 C have clustering structures. Although Morinaga [4] proposed the linear- 
chain structure for 12 C, the clustering states of 0^~(7.66 MeV) and 2^~(~10 MeV) have 
been discussed to have a loosely-coupled 3a configuration [19] rather than the linear- 
chain structure. We here investigate the stability of the 3a-linear-chain configuration 
in 12 C from a viewpoint of orthogonality to low-lying states by applying a microscopic 
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(b) 



Figure 1: The linear-chain configurations of 3a (a) and 4a (b). 

3a model. In this calculation, we use the Volkov No. 2 potential [18] with the Majorana 
value M = 0.592 as an effective nucleon-nucleon interaction. By using this interaction, 
the experimental binding energy of 12 C (7.4 MeV) from the 3a threshold, which is 
calculated as —83.69 MeV, is well reproduced. 

Various intrinsic 3a configurations can be described by a parameter set / 2 , 6), 
where li and l 2 are intervals between two a-clusters and 6 is a bending angle (See Fig. I 
a)). The linear-chain configuration corresponds to 6 = tt. Here, we treat li and 6 as 
generator coordinates and construct the basis states which have various sets of these 
coordinates. We employ li = 1,2,3,4,5 fm and 6 = tt ~ 7r/3. We regard l 2 as a 
variational parameter and apply a cooling method to search for the energy minimum. 
The 45 GCM basis states described by parameters (/ l5 / 2 , 9) are obtained, and their 
explicit values are presented in Table I. 

After projecting the obtained pure linear-chain intrinsic states (6 = tt) to J v = 
+ , we obtain the energy —71.72 MeV at the minimum point (/i, / 2 ) =(4 fm, 2.89 fm) 
(GCM basis No. 4). This energy is much higher than the 3a threshold where the 2 
state is observed. To solve the longitudinal vibration of three a clusters, we employ 5 
GCM basis states which are given by (/ l5 l 2 ) = (I fm, 3.55 fm), (2 fm, 3.22 fm), (3 fm, 
2.97 fm), (4 fm, 2.89 fm), and (5 fm, 2.92 fm), keeping 9 = tt (GCM basis of No.I ~ 
No. 5). After solving the GCM equation, we have the binding energy gain by 2.86 MeV 
( — 74.59 MeV), but its energy is still too short for the energy of the experimentary 
observed 0^ state. 

Furthermore, we take account of the bending motion (6 = ll7r/12, 107r/12 and 
97r/12) which corresponds to fluctuation around the linear-chain configuration. The 
solutions of the GCM equation with those basis states (No.I ~ No. 20) give the lowest 
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Table 1: The GCM basis states for 12 C described by l\ (fm), / 2 (fm) and 9 (radian) of 
Fig. 1 (a). / 2 is determined to minimize the energy at given li and 6. 

energy —80.40 MeV. Although this energy is about 3 MeV higher than the energy of 
the 0^ state, it seems to give a good correspondence with that. The wave function of 
this solution has a dominant components of the linear-chain configuration, which are 
seen from the fact that overlap with the linear-chain GCM solution (described in the 
previous paragraph) is 0.65. 

From the analyses described above, one may consider that the 0^ state can be 
interpreted as a 3a linear-chain state with small bending vibrations. However, this 
solution may also have a large overlap with the ground state of the 3a system, and it is 
necessary to solve the 0^ state by taking into account of the coupling with the ground 
state. In order to clarify this point, we performed the GCM calculation with larger 
basis states including those which describe low-lying levels. After diagonalizing the 
Hamiltonian with GCM basis states (No.l ~ No. 45 in table I), we obtain the ground 
state energy —91.23 MeV which well corresponds to the experimental one (—92.2 MeV). 

In Fig. 2, we show the conversion behavior of the energy spectrum. As the number 
of GCM basis states increases, each basis configuration generally becomes closer to 
equilateral triangle states. The ground + state is reproduced around the basis states 
No. 30 which corresponds to the configurations of 6 = 6tt/12. By this increase of 
basis states, the second + state converges energetically. This solution has the overlap 
with linear-chain (6 = tt) basis states of only 0.346 (absolute square), and cannot be 
interpreted linear-chain state. 
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Figure 2: The diagonalization of Hamiltonian of 12 C. 

3.2 4a states in 16 

In a similar way to the 3a system, we investigate the stability of the 4a linear-chain 
configuration in 16 0. We use the Volkov No. 2 interaction [18] with the Majorana pa- 
rameter M = 0.63, in order to reproduce the ground state energy of 16 ( — 127.6MeV) 
[17]. However, this effective nuclear interaction has a shortcoming that the calculated 
binding energy of 12 C (—82.0 MeV) does not reproduce the observed one (—92.2 MeV). 
This problem has been discussed [17] from both sides of model spaces and effective in- 
teractions. As one important conclusion, it has been shown that the clustering states 
are well reproduced around the corresponding threshold energy for various effective 
interactions. 

We apply the microscopic 4a-cluster model and carry out the GCM calculation 
for stability examination of the 4a linear-chain configuration. The employed GCM 
basis states are listed in Table II. The linear-chain states and their bending motions 
are described by 45 GCM basis states (No.l ~ No. 45 in Table II) of the 8 Be+ 8 Be 
configurations with length parameters (/ x and l 2 ) and an angle parameter (6 = tt ~ 7r/2) 
as shown in Fig.l (b). The relative position of two 8 Be clusters is determined by / 2 
{1 2 = 0.5,1.0,1.5,2.0 and 2.5 fm) and the distance (/i) between the two a-clusters of 
each 8 Be is chosen to minimize the energy at given / 2 and 6. In addition, we prepare 
15 GCM basis states obtained by the Constraint Cooling Method for the description 
of other lower states. 

As the first step in research of the 4a-linear-chain state, we search the + optimal 
linear-chain configuration with (/i = 3.29 fm, / 2 = 1.50 fm, 6 = 0) at the energy 
E = 34.79 MeV from the 4a threshold within a single Slater determinant. As the next 
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Table 2: The GCM basis states lor 16 0. Basis state configurations configuration No.l 
~ No. 45 are described by li (1m) / 2 (1m) and 6 (radian) ol Fig. 1 (b). The value ol / 2 
is determined to minimize the total energy. The GCM basis states No. 46 ~ No. 60 are 
obtained by the constraint cooling method, and classified with the constraint principal 
quantum numbers (n X} n y and n z ). 



step, we examine its stability lor mixing with configurations other than the linear-chain 
by solving the GCM equation with increasing the number ol the basis states. 

The + -state energy convergence is shown in Fig. 3, and the ground state energy 
converges to —129.1 MeV. The energy ol the + state obtained at —91.68 MeV (finally 
corresponds to Og ) converges with a very small number ol GCM basis states (= No. 15), 
which means this state contains large linear-chain components with small bending 
(< IOtt/12). As an evidence, this state has an overlap (absolute square) with the GCM 
basis states Nos.1-10 (6 = tt and Htt/12) ol 0.461 and with the GCM basis states 
Nos.1-15 (6 = tt, Htt/12 and IOtt/12) ol 0.613. Even when a larger number ol GCM 
basis states are adopted, the state ol —91.68 MeV is hardly affected by other lower 
states and approximately keeps its own energy value. However this energy is higher 
about 20 MeV compared with the 4a threshold energy ( — 111.6 MeV). 

There may be at least two reasons lor the too-high-energy solution ol Og : 1) The 
12 C+a threshold is not correctly reproduced due to the present NN interaction, and 
is calculated at the energy lower than experiment by about 10 MeV. As was discussed 
in the previous paper, it is very difficult problem to reproduce both threshold energies 
ol 12 C+a and 4a by an effective NN interaction. Because ol this shortcoming, 4a 
clusterized configurations are calculated at rather high excited energies. 2) The GCM 
basis states describing the Og may be insufficient. We have employed the configurations 
ol bending motion deviating Irom the linear-chain, and ol the three-dimensional 4a 
structure though they contribute little to this Og state. Then the energy ol this state 
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Figure 3: The diagonalization of Hamiltonian of 16 0. 

would be lower if we solve other small vibrations around the linear-chain configuration 
such as Z letter-type mode or three-dimensional twinst mode. 

However, an improved calculation based on the above reason for the fault would 
only increase the quantative precision, and at least we can suggest two things from the 
resent result. First, the state of — 91 .68 MeV has a large overlap with the linear-chain 
configuration which shows the 6 stability. Second, this state is hardly affected by other 
lower energy states and there is a large possibility to correspond to remain as resonance 
state. 

In Fig. 4, the linear-chain states with different angular momenta are presented. 
From the 2 + state, the rotational band is shown to appear with a very small slope of 
88 keV because of its large inertia moment. This moment of inertia is similar to the 
experimentally suggested value (about 60 keV). 

4 Summary and Conclusion 

In this peper, we have microscopically studied the existence and the stability of Net 
linear-chain state in 12 C and 16 0. We have emphasized the importance of the orthog- 
onality between the linear-chain state and other lower states. The linear-chain state 
of 12 C is hardly considered to be a stable state as we increase the number of the 
GCM basis states. In other words, their components fragment to the levels over broad 
range of excitation energy, especially to the Oj" state, and these couplings make a pure 
linear-chain state unstable. On the other hand, the linear-chain state of 16 is hardly 
affected by other states, and the rotational band structure is shown to appear, although 
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Figure 4: The rotational band of 16 linear-chain state. The numbers show the order 
of excited state. 

the absolute energy of the band head is higher by about 20 MeV compared with the 
4a-threshold. 

These conclusions are consistent with the previous work by Ikeda [11], where the 
stability against small vibrations around the equilibrium configuration is considered. 
However, we have shown in this work that the linear-chain state of 16 may survive 
even after the couplings to the lower states are taken into account. 
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